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We study the effect of varying the boundary condition on the spectral function of a finite ID 
Hubbard chain, which we compute using direct (Lanczos) diagonalization of the Hamiltonian. By 
direct comparison with the two-body response functions and with the exact solution of the Bethe 
Ansatz equations, we can identify both spinon and holon features in the spectra. At half-filling 
the spectra have the well-known structure of a low energy holon band and its shadow — which span 
the whole Brillouin zone — and a spinon band present for momenta less than the Fermi momentum. 
New features related to the twisted boundary condition are cusps in the spinon band. We show that 
the spectral building principle, adapted to account for both the finite system size and the twisted 
boundary condition, describes the spectra well in terms of single spinon and holon excitations. We 
argue that these finite-size effects are a signature of spin-charge separation and that their study 
should help establish the existence and nature of spin-charge separation in finite-size systems. 
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I. INTRODUCTION 

The nature of both the ground and excited states for 
the Hubbard model (HM) has been the theme of many 
recent studies. The HM itself is important not only be- 
cause of its simplicity but because it is believed to include 
many of the principle features of strong electron correla- 
tions in a wide range of real materials. These include the 
2D cuprate high-T c materials and in ID, compounds 
such as TTF-TCNQ [§ and SrCu0 2 @. 

The eigenvalues can be found exactly for the ID HM 
from the Bethe Ansatz (BA) equations 0. The corre- 
sponding eigenstates can be interpreted in terms of sep- 
arate populations of spin and charge excitations. The 
separate dynamics of spin and charge (spin-charge sepa- 
ration) together with other anomalous properties of the 
HM in ID Upl are now considered to be the standard 
signatures of non-Fermi (or Luttinger) Liquid behavior 
in ID systems and is found in other exactly soluble ID 
models like the Tomonoga-Luttinger model |7],|| . 

The exact BA solution is not available in higher- 
dimensional versions of the model nor even in ID when 
'extended' effects such as next nearest neighbor hopping 
are included. Instead, one has to resort to finite-size sys- 
tem calculations and approximate methods. One of the 
problems with finite-size system calculations is that the 
restriction to small systems can make information about 
the nature of low-energy excitations in the thermody- 
namic limit hard to extract. Even so exact diagonaliza- 
tions are an increasingly important method for studying 
strongly correlated electron systems || . 

Here we investigate how finite-size effects directly re- 
flect the nature of the underlying excitations. In partic- 
ular we look at how spin-charge separation manifests it- 
self in the variation of the spectral function as a function 
of the boundary condition. We work with the ID HM 



where the characterization of results of direct numerical 
diagonalization (in terms of spinons and holons) can be 
independently verified by comparing with the exact BA 
solution. 

There has already been considerable interest in how 
the HM behaves as the boundary conditions are changed 
by introducing a phase twist, S, in the periodic boundary 
conditions (PBC), ie c\ +N = e~ tS c\. This corresponds to 
threading the loop formed by the N sites in the chain 
with the magnetic flux <&o<5/27r, where $o is the flux 
quantum and the cj create electrons in orbitals at site 
i. One consequence of the phase twist is to shift all al- 
lowed momenta by an amount S/N yielding the allowed 
momenta: 



k = 



2mr + S 
N ' 



(1) 



for n integer. When an electron is added or removed, 
the many body system changes its total momentum by 
an amount corresponding to one of the allowed values in 
(Q). By varying S we can, in effect, simulate a continu- 
ously varying momentum in a finite system. 

The ground state energy of the ID HM has an inter- 
esting structure when plotted as a function of S [[ll] fl3| . 
For large U, there are oscillations in the ground state en- 
ergy which are often found to vary with more than one 
period. These oscillations — which may include periodic- 
ities of 1, 1/N e and N±/N e (for a total of N e electrons 
present, of them being spin-down) — give rise to cusps 
in the energy variation for a finite system. These effects 
are accounted for in terms of transitions mainly in the sea 
of spin particles as 5 is changed jll],[l2| . The ground state 
energy is found to be consistent with a hole in the spin 
sea which traverses from one Fermi point to the other as 
the boundary condition is varied. (Since the spin 'exci- 
tations' have zero energy, the spin sector can be used by 



1 



the system as a momentum bank charging no 'interest'.) 

To gain a quantitative picture of excitations in the ID 
HM, various authors have used different approximations 
and direct diagonalizations. By using Lanczos diagonal- 
ization, Kim and collaborators (|JT^] calculated the pho- 
toemission spectral function, A(k,u>), for a finite half- 
filled system with fixed PBC (they actually considered 
the t-J model). Their result, which is essentially the 
same as that of the strongly correlated HM, consists of 
spectral weight belonging to two bands, which they iden- 
tified with separated spin and charge degrees of freedom. 
The 'charge band' (of width 4\t\ and present in the whole 
Brillouin zone) traces out the dispersion of the so-called 
holon excitations. The 'spin band' (of width irJ/2 for 
J = 4t 2 /U) exists only for \k\ < kp (kp = ir/2) and is 
associated with an excitation of the corresponding spinon 
population. This structure of the spectral function is a 
natural consequence of spin-charge separation and can 
be explained through a semi-intuitive picture which has 
been termed the 'spectral building principle' (SBP) [jl5|j3| 
(see section IV below). 

In this paper we study the spectral functions of small 
ID systems (N — 12) at half-filling and in the presence of 
variable twisted PBC's. We show how the spectral build- 
ing principle can be adapted to account for the features 
that appear in the spinon bands and suggest that these 
features provide a natural test in finite system calcula- 
tions for spin-charge separation. The structure of this 
paper is as follows: In section II we outline the proce- 
dure which we adopt to calculate and analyze the spectral 
functions of the ID HM with twisted PBC. Our results 
are shown in section III. In section IV we review and 
adapt the SBP to help interpret our spectra. Finally in 
section V we summarize our findings and speculate how 
the SBP could have a bearing on systems other than the 
half-filled HM in ID. 



II. PROCEDURE 

We compute the spectral functions using the method 
used by Tohyama and Maekawa [jl6|. For a given 8 we 
compute the spectral function A(k,u), 

A{k,u) = 1 £ KV'm _1 |cfcil^*)| a *(w + [E^- 1 - <•]) 

m 

(2) 

at those momenta consistent with the choice of 8. We use 
the standard Lanczos recursion method p]]. In eq. (0), 
\4>m") ' s the mth energy eigenvector (with energy E^°) 
of the 7V e -electron system (at half-filling there are equal 
populations of each spin species) . All wavefunctions and 
energies are dependent upon 8. Following p6[ we ana- 
lyze the peaks in the spectral function by computing the 



complementary (fixed particle number) dynamic spin and 
charge correlation functions defined respectively as, 

% 1 w) = i£l(^ , " 1 l%)K , " 1 )l a 

m 

x^+K- 1 -<-!]), (3) 
N(q,uj) = ±J2M™~ 1 \ N ^0 e ~ 1 )\ 2 

m 

x<S( W + [£^-- 1 -<-- 1 ]). (4) 
The spin and charge operators, S(q) and N(q) are 

S (l) = ^2 exP^OK'T ~ n H ~ S o) ( 5 ) 

i 

N (l) = 5Zexp(igri)(« iT + n il - N ) (6) 

i 

where the summations are over sites, So an d No are the 
spin and charge backgrounds (So = (iVj — N±)/N and 
N = (7V T + Ni)/N) and n ia = c\ G c ia . 

The momentum transfer, q, is defined via q = kp — k. 
The Fermi momentum, kF is defined as the momentum 
state with the highest single electron energy which is oc- 
cupied at half-filling in the non-interacting (U = 0) case 
and which has positive momentum. The restriction to 
kp > is important in order to maintain the association 
of excitations around the positive Fermi point with small 
q. Consequently, for some choices of 8, {ipo"^ 1 ) m (Pf) 
and (0) is not the true ground state, but is the lowest en- 
ergy level in the subspace of states which share the same 
momentum as Ck F ilV^)- Since both k and kp depend 
upon 8 according to eq. (|l|), the effect of 8 cancels in q 
as expected {q is a momentum transfer and the spacing 
between the allowed wavenumbers k is fixed). 

The dynamic spin and charge correlation functions 
probe the properties of the excitations in the N e — 1 elec- 
tron system. A level with a strong response in S(q,Lu) is 
taken to indicate a spinon excitation. Features in N(q, u) 
are then associated with holon excitations. Some lev- 
els could be due to multiple excitations in one or both 
sets of particles. By comparing the spectral function to 
the corresponding spectra of S(q, ui) and N(q, to) we have 
been able to classify many spectral features by identify- 
ing which excitations are spinon-like and which are holon- 
like. This procedure is repeated over the range of 8- values 
(0 < 8 < 27r). This allows us to plot spectra as a func- 
tion of a quasi-continuous momentum variable, although 
strictly such spectra are taken from an ensemble of sys- 
tems with differing boundary conditions. 

Because the HM is integrable in ID, we can indepen- 
dently verify our assignments of the features we find in 
the spectral function by reference to solutions to the BA 
equations. For an A^-site Hubbard chain with ./V e parti- 
cles and 2Vj_ down spins and a twisted boundary condition 
these are (i = 1) 
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kjN = 2ttIj + 5-J2 2 



tan 



4 (sin ki — A Q ) 
U 



(7) 



and 



^2tan _1 



4(sin ki — X a ) 
U 



= 2irJ a 



(8) 



A', 



^tan 1 



2(A a - A/3 



The quantum numbers Ij, J a characterize the charge and 
spin degrees of freedom. They are either integer or half- 
odd integer according to the conditions 



/. = ^(modl) 

J a = z (modi). 



(9) 
(10) 



The total energy E and total momentum P of the system 
of particles is 



E = — 2 /J cos/cj 



P = y fc He 



6%. 

N 



(11) 
(12) 



The BA equations (Q) and (||) are easy to solve nu- 
merically assuming real rapidities kj , X a by simple root 
searching. 



III. RESULTS 



Figure [l] is the spectral function for the half- filled 12- 
site Hubbard band with U — 10 and antiperiodic (5 = 7r) 
boundary conditions computed using direct diagonaliza- 
tion (Lanczos). The curves are cross sections through the 
spectrum at those momenta consistent with this choice of 
5. The peaks have been Lorentzian broadened as a guide 
to the eye. In order to label the peaks (spin or charge- 
like), spectra of the spin and charge correlation functions 
have been superimposed and the dispersive tendencies of 
the main spin and charge-like peaks have been drawn in 
thin lines. The results are very similar to those of Kim 
et al. Ill 
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FIG. 1. The raw spectral functions computed using Lanc- 
zos for a ring of 12 sites at half-filling with U = 10 and 8 = n. 
The spectral functions (shaded) are shown together with the 
diagnostic spin (continuous line) and charge (broken line) cor- 
relation functions for the 11 electron system. The peaks have 
been Lorentzian broadened to Atu = 0.04|t| and the correla- 
tion functions have been scaled by a factor 5. Some peaks 
have been joined between different k values with thin lines to 
illustrate the dispersive nature of the principle excitations. 

To illustrate the essential features of the spectra ob- 
tained by exact numerical diagonalization as 5 is varied, 
we show the spectra in the k-u> plane. Each peak in 
the spectral function A(k,uj) is represented by a sym- 
bol centered about its momentum (k) and energy (lu). 
Figure ^a is for the same system as in Figure [j] but for 
a range of 5. The size of the symbols is used to de- 
note the spectral weight while (as indicated in the inset 
to Figure ||a) the shape is used to distinguish between 
spin and charge excitations. For some peaks denoted by 
shaded squares, there is a contribution from both types 
of correlation (as indicated for example in Figure [j] at 
k = 0.4177r). The remaining peaks which cannot be 
identified through the correlation functions are denoted 
by crosses. We also show the spin and charge corre- 
lation functions for this system (Figures |^b and ||c re- 
spectively) which were used to classify the features in 
the spectral function into 'charge- like' and 'spin- like'. 
(To make the comparison with A(k 1 uj), straightforward 
S(q = kp — k,uj) and N(q = fcj? — k,oj) are shown as 
functions of A;.) 

The vertical lines on Figure || denote those momenta 
corresponding to the choice of 8 in (Q) which give rise to a 
closed-shell configuration for the electrons in the ground 
state. For the 12 particle system this corresponds to 
5 = tt. The separation between two adjacent vertical 
lines corresponds to a change of 27r in 5. 
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FIG. 2. (a) The spectral function for the 12-site half-filled 
system with U = 10 computed by varying the boundary con- 
dition parameter 8. Peaks in the spectral function are clas- 
sified into 'spin-like', 'charge-like', 'both' and 'neither' using 
the spin and charge correlation functions for the 11 particle 
system shown in (b) and (c). The areas of the symbols is pro- 
portional to the spectral weight in (a) and to the amplitude 
of the correlation functions in (b) and (c). The continuous 
curves in 2(a) are those obtained from solving the BA equa- 
tions and the labeling a-f denotes the quantum numbers J a 
shown in Table M. 



In Figure 2, the solutions of the BA equations are su- 
perimposed on the results obtained by direct diagonaliza- 
tion. Each curve, labeled by the letters a-f, corresponds 
to a different choice of quantum numbers J a as shown in 
Table B. The Ij distribution in each case is that of a sin- 
gle holon, the position of which traces out the complete 
holon band as the phase shift is varied. The agreement 
between the labeling of the results of the direct diagonal- 
ization and the results of the BA calculation show that 



the numerical procedure based on comparison with spin 
and charge correlation is a very reliable way of identifying 
the nature of the excitations involved. 

In Figure |]a there are obvious features related to the 
twisted PBC. At low binding energies and low momenta 
(upper left part of Figure ||a), there is a series of 'humps' 
in the spin-like band, which have an approximate sym- 
metry about the closed shell lines. The main charge-like 
band, which exists over the full Brillouin zone (band 'a' 
in Figure |^a) , has the same shape as one would expect 
to recover in the thermodynamic limit. As we argue in 
section IV, these effects of the twisted PBC should be 
expected and are typical for systems with spin-charge 
separation. 

The extension of the charge band to negative momenta 
(not shown) and then reflected about k — gives the fea- 
ture which disperses with positive curvature in the high 
energy region of Figure ||a (band T). Something similar 
to this curve was seen by Penc et al. where it was 
interpreted as a shadow of the main holon band shifted 
in momentum by (we will refer to this band as the 
'shadow band'). Penc et al suggest that this shift is due 
to scattering by spin fluctuations at 2kp (the spin cor- 
relation function has a singularity at 2k p for strongly 
correlated electrons Jl0|]). Most of the shadow band is 
of undefined nature in our results (although some peaks 
do have a very weak spin correlation). The explanation 
of this feature in the context of the Spectral Building 
Principle (section IV) is compatible with that of Penc et 
al. 

There are some features in the photoemission spec- 
trum which our comparison with the two-particle cor- 
relation function does not clearly identify as being re- 
lated to either spin or charge correlations. These tend 
to have low spectral weight and/or to be at energies far 
from the Fermi energy. There are conversely also fea- 
tures in the correlation functions which are not relevant 
to photoemission. The spin correlation function for ex- 
ample exhibits very strong peaks at k > kp (Figure ^]b). 
These have no photoemission association except where 
they cross the charge band. We attribute these humps 
in the correlation function for k > kp to multiple holon 
and spinon excitations (section IV). 
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TABLE I. The spin quantum numbers for the various so- 
lutions (labeled a-f in Figures |^ and ^) to the Bethe Ansatz 
equations (^). The hole in the spinon band or 'antispinon' is 
denoted by X. 
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The new spinon related structure in the spectral func- 
tion is a finite-size effect. There are N/2 — 1 spinon cusps 
present for — n < k < ir. For this system and all systems 
with N/2 even, there is a cusp at k = 0, while for odd 
N/2 the cusp is displaced from k = by w/N. The range 
just below kp in which the spectral function has mixed 
spin and charge character is of width n/N and becomes 
a single point in the thermodynamic limit. 
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FIG. 3. As Figure |a but for U = 100 and (7 = 2. 

Results for U = 100 are shown in Figure^. The basic 
spectral structure of Figures §a and |pi is similar. There 
is more spectral weight in the high energy regions than 
for U = 10 with the shadow band having almost the same 
weight as the main holon band. The spinon band width 
parameter J ~ At 2 /U tends to zero in this limit and, in 
the infinite system, the spinon excitations for U — » oo 
do not disperse [pTsf| . Unlike the spinon bandwidth, the 
humps in the spinon part of the spectra do not scale with 
l/U. 

In the weakly correlated case of U — 2 shown in Figure 
, the main holon band for k > kp and the shadow band 
have a decreasing share of the spectral weight. (As the 
interaction is reduced, the spectrum is tending towards 
a conventional quasiparticle band with an electron Fermi 
surface.) There is still a two peak structure for k < kp, 



although the identification of spin and charge-like char- 
acter of the peaks in Figure ||b is only an indication of 
which excitation-type dominates, as there is notable mix- 
ing of spin and charge correlations. This is seen in both 
the charge and spin correlation functions which show a 
response at the corresponding frequencies. The humps 
in the weak holon band for k > kp develop into the one- 
particle two-hole excitations of the non-interacting case 
(where of course they have no spectral weight). 



IV. THE SPECTRAL BUILDING PRINCIPLE 
(SBP) 

The photoemission spectrum for a strongly correlated 
HM in the thermodynamic limit can be explained using 
the semi-phenomenological SBP The SBP deter- 

mines where spectral weight is expected and classifies the 
features as spin or charge-like. The SBP does not pre- 
dict the magnitude of the spectral weight itself and this 
information must be found by other means like exact di- 
agonalization or approximate methods. 

The SBP assumes that the ground state of the half- 
filled system is made up of a half-filled spinon band con- 
taining N± spinons and an empty holon band, with en- 
ergies E s (k s ) and Eh{kh) and momenta k s and kh for 
the spinon and holon respectively. The SBP treats the 
spinons and holons as independent excitations obeying 
the Pauli exclusion principle. At half-filling the disper- 
sion relations are given by |2C[| , 

E s {k s ) = |jcos(fc s ) (13) 
E h (k h ) = 2\t\cos(k h ). (14) 

The peaks on the spectral function are assumed to in- 
volve a single anti-spinon excitation (the annihilation of 
a spinon from the spinon sea) and a single holon excita- 
tion (the creation of a holon in the empty holon disper- 
sion curve). Taking the spinons and holons to be inde- 
pendent excitations, the electron Green's function would 
then be expected to be well-described by a convolution 
of the spinon and holon Green's functions. Conservation 
of energy and momentum would imply 

u = E s - E h (15) 
k = k s ~ k h . (16) 

Here u> and fc, as in Figures ||a and |L are the energy 
and momentum lost by the system when an electron is 
photocmitted. 

According to the SBP, a peak on the spinon band in the 
spectral function in the strongly correlated limit is inter- 
preted as a specific anti-spinon excitation in combination 
with the (fixed) lowest energy holon excitation at ftjj. In 
the thermodynamic limit, k® = it. By varying the choice 
of anti-spinon excitation and using (O)-(filf), this gives 



|b 



5 



the entire spinon band, w s (k) = (nJ/2) cos(fc + tt) + 2\t\ 
but with a cut-off at the spinon Fermi surface at k = 
±7r/2. In the photoemission spectrum, the occupied re- 
gions of the spinon dispersion have been shifted in energy 
and momentum by the holon excitation at fc°. Likewise, 
the SBP assumes that holon bands (main and shadow) 
are made up of all possible holon excitations together 
with a fixed anti-spinon excitation of lowest energy at 
kg = ±7r/2, corresponding to spinon annihilation at the 
spinon Fermi points. The holon bands in the spectral 
function have the form, u>h(k) = —2\t\ cos(fc=F7r/2). Since 
the holon bands measure the unoccupied regions of the 
holon dispersion, the holon band observed in the photoe- 
mission spectrum extends throughout the Brillouin zone. 

The SBP as outlined in [|l5]]|] describes the spectra in 
the thermodynamic limit. In order to proceed to finite 
systems and, in particular, to explain the effect of the 
twisted PBC we need to determine how the spinons and 
holons should respond to 8. We assume that the holon 
and spinon momenta are given by 



2?r7,- 



k - 



N 
2-kJq. 



N 



(17) 
(18) 



with the Ij and J a integer or half-integer as in (g) and 
(|io|). These quantities were introduced previously in ||] 
to help elucidate the connection between the exact eigen- 
states of the Hubbard model and the exact solution of the 
Tomonaga-Luttinger liquid. We have generalized them 
by including the phase shift term in (|i~7|). 

The allowed values for the holon quantum numbers 
kh, k s or (equivalently) Ij, J a are functions of the num- 
ber of particles N e and the spin polarization iVj_. The 
photoemission process measures the energy of the N e — 1 
particle system as a function of the missing energy and 
momenta, so it is natural to expect the quantum numbers 
of the N e — 1 particle system to characterize this miss- 
ing energy and momentum. We therefore assume that 
the parity conditions on Ij and J a are those for the the 
N e — 1 particle system. 

When we compute the spectral function using Lanc- 
zos and the energy dispersion by directly solving the 
BA equations, we measure frequency with respect to the 
ground state of the system for the same boundary condi- 
tion. This energy itself depends on 8 and this effect is not 
accounted for by the SBP when we assume the relation 
(|l5|). For U — 2, the variation in the ground state en- 
ergy of the half-filled case is less than 10% of the spinon 
bandwidth and is smaller for larger U. For the spectral 
function taken for systems which are not half-filled and 
with no gap the picture is more complicated and will be 
addressed in a later publication. 

Figure |] shows the positions of the spinon and holon 
peaks computed using the SBP for the case U = 100 and 



should be compared with the numerical results of Fig- 
ure |a. The spinon regions are denoted by squares and 
the holons with circles. The agreement is good bearing 
in mind that only the spectral features which are asso- 
ciated with and are reproduced in Figure |^) . The 
main holon band, the shadow band and the spinon cusps 
are reproduced. The number of cusps and their positions 
agree well with the numerical result and this picture also 
gives an offset of the top-most energy of the holon band 
to a slightly lower momentum than fcp which is evident 
in the numerical result too. Shown for comparison is the 
result for the thermodynamic limit where there are no 
cusps and the top of the holon band is at fc^ exactly. 



(7=100 




FIG. 4. The main picture (left) shows the photoemis- 
sion spectrum for N = 12, N e = 12, N± = 6 and U = 100 
computed on the basis of the SBP (squares=spinon, cir- 
cles=holon) and should be compared with the result of Figure 
gp. The shaded points are for reference in the text to follow. 
The small plot shows the corresponding result in the thermo- 
dynamic limit (continuous line=charge, broken line=spin). 

The shape of the spectra Figure [| can be understood by 
considering the holon and spinon dispersions for a given 
8 shown in Figure For the spinon band in Figure ^ 
the set of occupied states of the spinon dispersion are all 
shifted together by the lowest energy holon excitation — 
labeled 'A' in Figure p| (note that the spinon bandwidth is 
exaggerated in Figure^). According to (13) and ( fL4[ ) this 
gives one set of equally spaced points on the spinon band 
in Figure [| (the shaded points) . Varying 8 changes both 
the energy and momentum of the lowest energy holon 
while the spinon energy is (by assumption) independent 
of 8. As a result, each peak in the spinon band of the 
spectral function is shifted to a slightly different position. 
The humps observed in the spinon band of the spectral 
function then simply reflect the variation with 8 of the 
holon state labeled by 'A'. At certain values of 8 the low- 
est energy holon excitation changes its momentum dis- 
continuously as the state previously labeled by 'A' moves 
up the holon band. Whenever this occurs it gives rise to 
the cusps in the spinon band in the spectral function. 



6 



HOLON SPINON 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

k h ln k s l 7i 



FIG. 5. The holon and spinon dispersions for the fi- 
nite system N = 12, iV e = 12, N L = 6 and 8 = 1.6tt. 
Open (filled) symbols denote unoccupied (occupied) holon 
and spinon states. The allowed holon states change with 5 
but those of the spinon are static. 

The form of holon bands in the spectral function can 
also be understood using the SBP. The shift of the holon 
dispersion by the spinon excitation 'B' gives rise to the 
main holon band (the upper holon band in Figure 
Since (by assumption) the allowed momenta in the spinon 
dispersion are independent of 5, the holon dispersion is 
shifted by a constant amount as S is varied. This is the 
reason why the holon bands in the numerical spectra are 
smooth and free of cusps in the strongly correlated limit. 
The alternative spinon excitation, 'C, has the same en- 
ergy as 'B' but shifts the holon dispersion by a different 
momentum. This leads to the so-called shadow band in 
the spectral function (lower holon band in Figure ^) . The 
momentum difference between 'B' and 'C is 2kF (true 
for all fillings) which explains why the shadow band is 
displaced by 2&f from the main holon band. 

The coupling between spin and charge which exists for 
moderate U means that assuming completely indepen- 
dent spinon and holon bands as in @ and is no 
longer strictly valid. Even so it is still possible to ex- 
plain the spectra using the SBP, although the non-zero 
spinon bandwidth J introduces a minor complication as 
we explain below. 

If we implement the SBP exactly as described above 
for the large U case, the low energy peaks predicted by 
the SBP do not agree completely with what we find us- 
ing direct diagonalization and the BA solution. Figure 
shows the photoemission spectrum computed on the 
basis of the SBP for U = 10 (J = 0.4) and N = 12 
at half-filling (shown in the inset is the N — > oo result 
for comparison). Unlike the numerical result of Figure 
there appear to be discontinuities the spinon band 
while the holon bands look very similar to what we found 
numerically (see Figure ^a) . These discontinuities come 
about because we have assumed that, as in the large U 
case, the spinon band of the photoemission spectrum cor- 
responds to the combination of a spinon with the lowest 
energy holon excitations as before. Whereas this is cor- 



rect for most momenta (with k < n/2) in Figure ga, it 
is not true for the range of momenta highlighted by the 
horizontal 'bars' in the top-left of Figure ||a. We find that 
for these electron momenta, lower energy electronic states 
can be formed by combining a spinon excitation with the 
slightly higher energy holon excitations at point 'D' in 
Figure g. Figure gb shows the photoemission spectrum 
computed using the SBP adjusted to allow also for low 
energy spinon and holon combinations (ie using states at 
'D' as well as 'A' in Figure |5|). Although it might seem 
peculiar that an electronic state formed from a higher 
energy holon state ('D') can be responsible for lower en- 
ergy electronic states, it can be explained as follows. The 
holon at point 'D' in Figured i^h> ^h) nas a higher mo- 
mentum than 'A' (kfi,Efi) (after adding 27r to fcf ). To 
form a hole with given momentum, the spinon state com- 
bined with the holon state at ! D' will be at a different 
point on the spinon band. The cost in energy of using the 
more energetic holon state 'D' can then be compensated 
by a lower spinon energy. This adjustment is not neces- 
sary for the U — ► oo case since there the spinon band is 
flat. 



(a) !/=10.0 (b) (7=10.0 
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FIG. 6. Spectral function computed using the spectral 
building principle for N = 12, N e = 12, N l = 6 and U = 10. 
In (a) the spinon (holon) band is shifted by the very lowest 
energy holon (anti-spinon) excitations only. The inset is for 
the thermodynamic limit (continuous line = charge, broken 
line = spin). In (b) we allow the spinon dispersion to be 
shifted by the next lowest energy holon excitations also. This 
should be compared with Figure ^a. 

Although U = 2 does not represent a weakly correlated 
system, it is too small for the SBP to work well. We find, 
in fact, that the SBP ceases to work well below U w 4, 
where the holon band starts to develop humps as a func- 
tion of 8. We have tracked these humps in the exact BA 
solution as a function of U down to small U (U < 0.1). 
They develop into the two-hole and one-particle excita- 
tions of the non-interacting system, where of course they 
have vanishing spectral weight. 

In the SBP plots (Figures 

| and g), we have shown 
only those peaks which are firmly classified using the re- 
sults of the spin and charge correlation functions. These 
turn out to those associated with the lowest (or in some 
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cases low-lying) spinon or holon excitations. The nu- 
merical spectra contain many other peaks in addition to 
the spinon and holon bands. By taking every possible 
combination of single holon and anti-spinon excitations - 
regardless of energy - the SBP can reproduce essentially 
all the peaks in the spectral function. It is clear then 
that the spectral function are dominated by the creation 
of single holon plus anti-spinon excitations. 

Other classes of excitation exist which can be ac- 
counted for using the SBP but which are invisible in the 
spectral function. For example for the U — 10 system 
of Figure it is clear that that there are contributions 
to the spin correlation function for k > kp for momenta 
k > 7r/2 (there is no significant charge correlation in 
addition to that already discussed). We attribute these 
'extra' humps to multiple excitations, namely the holon 
and anti-spinon excitations already described accompa- 
nied by a spinon/anti-spinon pair (a spin wave). The 
effect of the twisted PBC makes these easy to identify, as 
the spinon humps at k > kp appear as images of one or 
more of the original humps in the spinon band but shifted 
by the extra momentum and energy of the spinon/anti- 
spinon pair. The spin waves can provide potentially any 
momentum in the set {2nit /N} and an energy which is 
finite but small (due to the narrow width of the spinon 
band) . 



V. SUMMARY 



We have studied the photoemission spectrum of ID 
Hubbard model at half- filling as a function of the periodic 
boundary condition (PBC) in finite size chains. We have 
shown that the finite-size effects can be well explained us- 
ing the spectral building principle (SBP, j||lfj) adapted 
to the case of finite systems and the twisted PBC, in 
terms of independent populations of spinons and holons. 
We have found that the photoemission spectral functions 
are explained for U >^ 4 by assuming that only the 
holons are directly dependent on the twisted PBC and 
that one holon-one antispinon pairs account for essen- 
tially all the spectral weight. Other excitations involving 
additional spin-wave processes are observed in the spin 
correlation function. 

The key feature found in the strongly correlated sys- 
tems (U >^ 4), in which the spinons and holons may be 
thought of as independent excitations, are the 'humps' in 
the spinon bands of the photoemission spectra. We have 
shown that the curvature of these humps directly reflects 
the curvature of the holon dispersion at the bottom of 
the holon band and so could be used to extract directly 



the holon effective mass in cases where there is no direct 
independent way to compute this quantity. 

We have restricted our study to the ID nearest neigh- 
bor Hubbard model where we have been able to show 
that identifying the spin and charge excitation bands us- 
ing the corresponding correlation functions agrees exactly 
with the results from Bethe Ansatz exact calculations. 
These results suggest that this method can be extended 
to other ID systems and higher dimensional systems for 



which there is no exact BA solution for the model |21|. 

We wish to thank the EPSRC for providing the funding 
for this work. 
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